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Abstract 

This papei0 gives a geometric description of functional spaces re- 
lated to Domain Decomposition techniques for computing solutions 
of Laplace and Helmholtz equations. Understanding the geometric 
structure of these spaces leads to algorithms for solving the equations. 
It leads also to a new interpretation of classical algorithms, enhanc- 
ing convergence. The algorithms are given and convergence is proved. 
This is done by building tools enabling geometric interpretations of the 
operators related to Domain Decomposition technique. The Despres 
operators, expressing conservation of energy for Helmholtz equation, 
are defined on the fictitious boundary and their spectral properties 
proved. It turns to be the key for proving convergence of the given 
algorithm for Helmholtz equation in a non-dissipating cavity. 
Using these tools, one can prove that the Domain Decomposition set- 
ting for the Helmholtz equation leads to an ill-posed problem. Never- 
theless, one can prove that if a solution exists, it is unique. And that 
the algorithm do converge to the solution. 

1 Introduction 

In the framework of domain decomposition, given a bounded open set Q = 
Hi U 02 U T where the two open sets fli and Q2 are not overlapping, and T a 
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common subset of their boundary (called the fictitious boundary), and given 
a (global) solution u of the Helmholtz equation 

Au + k 2 u = f e L 2 (n) and u e if* (ft) 

the aim of this paper is to understand the dynamics of the sequence (v ™, v J ) ne jv 
solving separatly the Helmholtz equations on fti and f2 2 , when equating the 
fluxes through T: (to — 1, 2 resp. to/ = 2, 1) 

° V ra _ ■ n _ UV m< _ ■ n-1 nn p 

Its ultimate aim is to prove convergence to (u^, u\n 2 ) of the sequence (u™, u^neJV 
solving the Helmholtz equations on Q\ and with a penalization on the 
boundary T is added, namely: 

For this sake, the geometry of the set of solutions of the Helmholtz equation 
on Q,i x Q 2 with equated energy fluxes is studied, through the study of the 
coupling operator defined on £ 2 (r) x £ 2 (r) which intertwins the fluxes. It 
turns out that the key for understanding the convergence of the sequence 
(ui,u%) ne N is the analysis of the spectral properties of the intertwinnig op- 
erator. 

Using these tools, one can prove that the Domain Decomposition setting 
for the Helmholtz equation leads to an ill-posed problem. Nevertheless, one 
can prove that if a solution exists, it is unique. And that the algorithm do 
converge to the solution. 

Convergence of the penalized algorithm is proven and numerical tests for 
solving the Helmholtz equation through this domain decomposition algorithm 
are given. 

The geometric analysis given here provides the theoretical background for 
another numerical algorithm for computing the global solution u, by a specific 
spectral method. A forthcoming paper describes and gives the numerical 
analysis of this algorithm. 

This domain decomposition algorithm (in a dissipating cavity case, i.e. with 
a Sommerfeld-like radiation condition on part of the boundary), was first 
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initiated and studied by B.Despres in [Dl] [D2] [BD], and computational 
results given by J.D.Benamou [B] [BD], F.Collino and P.Joly [CGJ]. 

In order to perform the geometric analysis of the set of solutions of the 
Helmholtz equation on QixQ 2 , one has first to make a complete description 
of the geometry of the set of solutions of the Laplace equation on fii x Q 2 - 
Geometric properties of this set proven below makes it possible to revisit 
the classical penalized Dirichlet/Neumann domain decomposition algorithm 
(with penalization) for solving the Laplace equation. A new version of this 
algorithm is given here, and proved to converge to the global solution, en- 
hancing the usual assumption on the penalization parameter. 

This completes classical results by O.Widlund [PW], P.L. Lions [L], or A.Quar- 
teroni and A.Valli [ FMQT] [FQZ] [QV]. 

The paper is organized as follows: in section 2 basic facts are revisited, 
although classical, and completed in order to set the geometric framework 
needed. (It also makes the paper self contained). A precise study of duality, 
and the link with the Poincare-Steklov operators, is performed, which turns 
to be central for the remainder of the paper. In section 3 a new version of 
the Dirichlet/Neumann algorithm for the Laplace equation is given, and con- 
vergence is proved. In section 4 geometric tools for the Helmholtz equation, 
and related domain decomposition algorithm, are given. Despres operators 
are studied and their spectral properties investigated. As is the intertwinnig 
operator. In section 5, convergence of the domain decomposition algorithm 
for Helmholtz equation is proved. In section 6 numerical tests are given. 

Throughout this paper, when dealing with the Helmholtz equation, the fre- 
quency k is assumed to be non-resonnant for the Dirichlet boundary condi- 
tion. More precisely we shall always make the following 
Assumption (A) —k 2 is not an eigenvalue of the Laplace operator on f2 
with Dirichlet boundary condition, i.e. the following problem is well posed 
for feL 2 (n): 

Am + k 2 u = f and u e Hq(Q) 
We shall also adopt the following 

Notation (N) normal derivatives at the boundary of an open set are always 
meant as the derivative along the outward unit normal vector 
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2 Basics 



Let Q C R d be a bounded open set whose boundary dVt is a C 1 -submanifold 
of _R d . Let T be an open C°°-submanifold of R d , such that: 

i] = Q 1 ufi 2 ur, orii = (dn n u r, <9fi 2 = (9ft n dn 2 ) u r 

where Qi and are open sets in R d . We assume that Q± and fulfill the 
strict cone property (see [Ag] for instance) and that T is transverse to dfl in 
the following sense: T is a C^-submanifold of R d with boundary, and there 
exists a < 1 such that for any a G dil fl T, we have: 

- a < n T (a).n d n(a) < a (1) 

where np(a) G C°(r) is a unit vector normal to T at a and ng^(a) G C°(dQ) 
a unit vector normal to dQ at o". 



2.1 Functional spaces associated to T 

Let Hq(Q) be endowed with the scalar product 



[u,v) H i {n) = J^VuVvdx 



For m — 1,2 let if m = {w G iJ 1 (r2 m ); ii|an n an m = 0}- Boundedness of the 

trace operators from ifg (Q) to H§ (dVt fl imply that these are Hilbert 

spaces when endowed with the scalar products: 

(u,v) Hm = VuVvdx 



Let p r (resp. p v m for m = 1,2) be the trace operator on T, i.e. the bounded 
linear operator from Hq(£1) (resp. H m ) to H l l 2 (T) which maps u to u\r- 
Let 

A = {u {r ;ue H^(Q)} = Hl{tt)/Kerp T ~ (fsTerp r ) ± 
and for m = 1, 2 

A m = {u\ T ; u G # m } = H m /Kerp T m ~ (Xerp^)- 1 
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Remark 1 Obviously Kerp T m = H%(n m ), A m C H l l 2 {T), A C H 1/2 (T) 

Because p r and p v m are bounded, A and A m are Hilbert spaces when endowed 
with the following norms: 

VA G A, ||A|| A = inf ||«||jfi(n) and VA G A m , ||A|| Am = inf \\u\\ H „ 

{u;ri|r=A} 1 {tt;ii|p=A} 

Remark 2 Obviously, for any v G iJg (^) and w G 

II/OOIU < \M H i {Q ) and ||p^H|| A < \\w\\ Hm 

Proposition 1 Let m = 1,2. For any A G A m , /x G A m : 

1- There exists a unique G H m such that Au^ = in Q m and 

pL«) = A 

0- One has: \\X\\ Am = \\u^\\ Hm and (A,/i) Am = (u^,u^) Hm 
3- For any A G A let u x = on Q m , m — 1, 2. T/jen 

||\||2 _ II A||2 , || A||2 
ll A llA — \\ u l lli/j + ll M 2 HfTa 

^- Using the previous notation, for any A G A, /i G A 

(A,/i) A = (M A ,M M ) H i(n) 

proo/: 

1- Uniqueness follows well posedness of the Laplace problem in i?g(fl m ). 
In order to prove existence, let u G be such that A = p^iu). Then 
Au G if _1 (f2 m ). Let v be the unique solution in H^Clm) of At> = Am. Then 
u^ n = u — v fulfills the property. 

2- Because of remark [TJ one has to show: 

which follows from the Green formula. 

3- Because p\(u x ) = p\{u x ) = A one has u x G Hq{VL). Obviously it is 
orthogonal to Kerp r . so 

H \ ||2 _ II A||2 _ || A||2 , || A II 2 

II a IIa — \\ u WmiQ) — \\ u i\\h 1 + \\ u 2\\h 2 
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2.2 A = Ai = A 2 

By symmetry, it is enough to prove A = Ai. In order to prove this alge- 
braic and topological equality, two key tools are needed. The first tool is 
the Calderon extension theorem [Ag], which applies here because Qi has 
the strict cone property, by assumption, and which gives a bounded linear 
operator E from if 1 (fix) to H l (R d ) such that: 

G H 1 ^), Ew {ni = w 

The second key tool is: 

Theorem 1 There exists a bounded linear operator r in H 1 (R d ) such that: 
\fv G H^B*), v mindn = {rv\ dn = and ru| ni = v\ ai ) 

proof: 

Assumption ([T]) gives a finite open covering {ujj)j of T n dVL and a change of 
variables (a?)j such that Vj = a?(u)j) is a neighbourhood of zero in R d and: 

a^THujj) = {z j G Vj] 4 = 0,4> 0} and ^(dQnuj) = {z j G Vj] 4 = 0} 

a j {Qi n ujj) = {z j G Vf, z{ <0,z J 2 > 0} 
a j {Q 2 n a,) = {z j G Vf, 4 > o, 4 > °} 

Regularity of the submanifold T gives an open covering {oj' k )k of Y and change 
of variables {b k ) k such that Wk = b k {uj' k ) is a neighbourhood of zero in R d 
and: 

b k {TnJ k ) = {z k G W k \z\ = 0} 

n u^) = {z k G W fc ; 2* < 0} and b k (Q 2 n w£) = G W & ; > 0} 

Compactness of T enables to select a finite sub covering of T still denoted by 
(ujj)j U {oj' k )k having the previous properties. 

Let to 1 = R d \ Hi and to 2 = R d \ tt 2 , so: 

R d = u 1 U u 2 U (UjWj) U (U fc c4) 

Let 

(a 1 , a 2 , (ay),-, (a' k ) k ) 
be a C°°-partition of unity associated with this open covering of R d . 
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Let e > be such that \/z G (UjVj) U (U k W k ) < z{ < e ^> z E tt 2 

Let ip(s) G C°°(R) be equal to one for s < and zero for s > e. 

Let ip{6) G C°°(R) be equal to zero for 6 < and equal to one for 6 > | 

For any u G H\R d ), 



we define tv as: 



r(« 1 f ) + r(a 2 v ) + T^jf) + X! r ( Q i' y ) 



with: 

T^V) = 

r(o; 2 t>) = a 2 v 



r(a' k v)(z k ) = ifj(z k )a' k (z k )v(z k ) in the local coordinates. 
These three quantities are multiplication of v by C°° functions, which are 
bounded as well as all their derivatives. It is linear and bounded in H 1 (R d ) 
with respect to v G H l (R d ). 

In order to define r(aijv), we first write otjV in the cylindrical coordinates as 
follows: 

ajv(z{, z 3 2 i z J 3 , .., z d ) = Ojvij^ , 6\ z%, .., z 3 d ) with z{ = r- 7 cos(^ J ), ^ = r^sin{6^) 
and define T(a:ji>) in these coordinates as: 



This quantity is linear with respect to v, and we prove its boundedness in 
H l (R d ) with respect to v G as follows (we omit the index j and 

denote the measure dz s ...dzd by dz): 



r( aj v)(r\ P, zl .., z d ) = ip{e j )&rv{ri , P, 4, .., z d ) 



t(ocv)\\ 



2 



r(av)\ 2 rdrd9dz + 



f | ^ 



■r(av)\ 2 rdrd9dz 
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+ / — \ip(9)— av + (p'(9)cn)\ 2 rdrd9dz + / |<£(0)-^ — av\ 2 rdrd9dz 
Jv r l 69 3 7y ctej 

<swp|(^| 2 [/ \cw\ 2 rdrd9dz + [ \-^-av\ 2 rdrd9dz 
Jv Jv or 

+ / ^\^6tv\ 2 rdrd9d~z + / |7^5w| 2 r<ir<i#od;l 
r 2 09 3 azj 

+2sup|<//| 2 / — \3v\ 2 rdrdOdz 
Jv r 2 

< 2sup\(p\ 2 \\av\\ 2 H i( Rd -) + 2sup\ip'\ 2 J — \av\ 2 rdrd9dz 

In order to estimate this last quantity we use the assumption v\qq iU qq = 
to have: 



av(r,9,z) = — —av(r, s, z)ds 
Je 09 

which gives (with e' the radius of the support of a in the r variable, and B 
a ball containing the support of a in the z variable): 



r I r f e ' p 71 1 p n 

/ — \av\ rdrd9dz < / / / — I / —av(r, s,z)ds\ 2 rdrd9dz 
Jv r 2 Jb Jo J-tt r 2 Je 09 

r re' rir j Q 

< 47r 2 III — \—c^(r,s,z)\ 2 rdrdsdz<4n 2 \\av\\ 2 H i (Rd) 
Jb Jo J-tt r 2 09 y ' ' 

we summarize to have: 

\\ T i.^)\\ 2 m{Rd) < C\\av\\ 2 HHRd) < C'\\v\\ 2 m{Rd) 

and this ends the proof of the boundedness of r in if 1 (_R d ). 

We end the proof of theorem [T] using the following obvious observations: 
T v \cii = because ip = 1 on R_ and if = 1 for 9 > | 
™|9Qiu9c = by assumption 

^idf^uac — because cp(0) = and ip = for z{ > e. 

Corollary 1 For m = 1,2, there exists a bounded linear map E m from H„ 
to Hq(Q) such that 

Vw G H rn (E m u)\ Qm = u 
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proof: for m = 1 for instance let E\u be the restriction to Vt of T-E/u.Boundedness 
follows from theorem [TJ 

Corollary 2 Ai = A 2 = A and the three norms ||.||a, H-IUx and ||.||a 2 are 
equivalent. 

proof: obviously A C A m and the previous corollary gives the converse inclu- 
sion. Moreover we have ||.||a > IUU m and the previous corollary gives the 
converse inequality. 

Corollary 3 T>(T) is a dense subspace of A for any of the three norms. 

proof: T>(T) is the set of traces on T of functions in V(Q) because T is a 
C°°submanifold. Density of T>{Y) for the ||.||a norm follows density of V(Q) 
in Hq(Q). Equivalence of the three norms ends the proof. 

i i 
Remark 3 If we denote as usual by Hq (T) the closure of V{Vt) in Hz(T) 

( which exists because V is C°°), then the previous corollary asserts that A C 
i 

Hq(T). Boundedness of the trace operators gives constants C,C m such that: 
VA G A, IIAII i < CIIAIIa and VA G A m , , IIAII i < C m ||A|| A 

2.3 Well-posedness of the Laplace-Dirichlet problem 

in H~ 1 (Q m ) x A 

Theorem 2 For m = 1, 2 

1- for any (f, A) G H~ l {yt m ) x A there exists a unique u G H\ such that 

Au = f in VL m and Pm( u ) = ^ 

2- we have the estimate \\u\\ Hm < \\f\\ H -^n m ) + ||A|| Am 

proof: Let be given by Proposition [TJ Let v = u — u^. The problem is 
equivalent to 

v G Hl(VL m ) and Av = f 

This is a well-posed problem and we have, because the Riesz representation 
operator is isometric, \\v \\H^(n m ) — WfWH- 1 ^)- 

So 

Nlflm < \Wm\\H m + \\v\\ Hm < \\f\\H-i(n m ) + ||A|| Am 
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2.4 Duality and the Poincare-Steklov operators 

Let A' denote the dual space to A, endowed with one of the three equivalent 
norms associated with the equivalent norms on A defined previously. 
We denote by (., .)aa' the duality product, and by (., .)vv the duality product 
in V'(T). 

Because of Corollary [3] and remark [31 we have the usual injections: 
PcAc hI(t) c L 2 (T) c H~*(r) cA'cV 

and for any A G T>(T) and v G A': (v, A)aa' = { u i X)vv> = [y, A)l 2 (T) 
where rj denotes the complex conjugate of functions or distributions r\. 
For any of the three scalar products on A we have (A, 77) = (A, rj) so for all 
norms: 

VAgA, ||A|| = ||A|| and W G A', \\v\\ = \\v\\ 

Notation 1 Let S denote the antilinear Riesz representation operator for 
A' in the A-scalar product, and S m (m = 1,2) this representation in the 
A m - scalar product, i.e. 

VA G A, Vz/ G A' (u, A)aa' = (A, S~ V) A = (A, S^u) Al = (A, S^v) M 

Let S (resp. S m ) be the linear isometric bisections from A to A' defined by 

VAgA SX = SX, S 1 X = S 1 J, S 2 X = S 2 J 

We denote by n m the normal unit vector on T pointing outward with respect 
to Q m . We denote by -J^f- the normal derivative on T of (p G T>(Q,), and by 
-J*— bounded extensions of this operator to any functional space. 

We use Proposition [1] to have: 



Proposition 2 



VA G A S m \ = and S = S 1 + S 2 

on m 



proof: by Green formula 

VA G A, n G A [S m r], A) Am A' m = (A, r,) i 
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Jil m JF (J I I'm 



rn 

I 



this proves that the distribution is bounded in the A m norm, so |^ G A 
and 

~ _ dull _ du m _ du m 

JmV "o "o "o 

dn m dn m dn m 

Corollary 4 

jS'A = SX S m X = *S* m A 
Remark 4 For any n G if m sncn toai An G L 2 (fi m ), we nane G A' and 

<9n 

||^-||A'<C(||n|| ffm +||An|| L2(nm) ) 

This is because for any ip G V(fl) we have: 

, dv 



■ ■> <4>\t)w = I V(pVvdx+ / (fAvdx 

•J Q jn J 



and this formula shows that the distribution is bounded on A. 

2.5 Adjoints 

Notation 2 : 

1- For a bounded linear operator T from A to A', we denote by T' its adjoint 
for the (A, A') duality, i.e. 

VA G A, Vn G A (Tn, A)aa' = (T'X, 77) aa' 

2- For a bounded linear operator T from A to A ; we denote by T* the adjoint 
operator in A, i.e. 

VA G A, Vn G A (Tn,A) A = (n,T*A) A 
Proposition 3 Form = 1,2 

S' — S S' m = S m 
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proof: by definition of S we have: VA 6 A, Vt/ G A 

(SX,rj) AA > = (SX,r]) AA , = (r?, A) a = (Kv) A = (SrJ,X) AAI = ( 5r ?' A )a,A' 

Theorem 3 For m = 1, 2 Zet m 7 = 2,1. For all A G A, 7/ G A 

1- (S^-Sm'V, A) A™ = (77, A) Am , 

2- (S~ S m 'T], A) Am , = (5 m /T7, 5 , m /A) A ^ 

5- S^S m > is self adjoint in A m and in A m i 
4- (ST 1 ^ + Si) is selfadjoint in A 

proof: 

1- (S^-Sm'V, A) Am = (A, S- l S m ,r}) Am = (S m ,rJ, A) AA ' = (A, rj) Am , = («, A) Am , 

2- We use Proposition [3] to have 

(S^Sm'Vi tyA m , = (Sm'X, 5'~ 1 S' m /?7) AA ^= (Sm'T/, 5'~ 1 5' m /A) AA / = 

{SmSmSm'i], S'~ 1 S' m /A) AA / = (5'~ 1 5' m /A, S^-S m iTj) Am = (S m <r], S m 'X) A ' m 

3- We use Proposition [3] to have 

(S^Sm'Vi A) Am = {S m 'fj, A)aa' = (^ m 'A,r7) AA / = (rjiS^Sm'X^Am = (v, S^S^X)^ 
On the other hand 

(S^Sm'T], A) Am , = (S^'A, S^Sm'T])^ = (S m 'T], S^ 1 S m ' X) AA > = (S~ S m > A, 7/) Am , = 

4- follows 3 

Corollary 5 Coerciveness: For m = 1,2 let m 1 = 2,1. There exists C > 
s«c/i i/jcrf /or a// A G A: 

1- (^- 1 ^A,A) Am = ||A||i m/ >C||A||i m 

2- {S m 'S ml X,X) Aml = \\S ml X\\\, m >C\\X\\l m , 

3- ((S{ 1 S 2 + S 2 1 S 1 )X,X) A > (l + C)||A||i 

2.6 On the Neumann problem 

Proposition 4 For m = 1, 2 and for any v G A' £/iere exists a unique u m G 

Au m = in fi m an<i m = v on T 

dn m 

Moreover 

\\ u m\\H m = \W\\a' 



12 



proof: uniqueness is straightforward, and existence is provided by proposition 
[T] and u m = u Sm v 

Moreover isometry of the Riesz representation gives: 

\\u m \\ Hm = KMU = \W\\a> 

Proposition 5 For m = 1, 2 and for any f G L 2 (Q m ) there exists a unique 
u m G H m such that 

du 

Au m = f and — — — =0 on F 

Moreover 

\\u m \\H m < C||/IU2 ( n m ) 

proof: Uniqueness is straightforward. For existence let v m G HQ(Q m ) be the 
unique solution of Av m = f. In remark H] we have shown that v = G A' 
and 

Hk < C(||/|U 2( n m) +||^|U m ) < C(\\f\\L*(n m) +\\f\\ H -Hn m )) < C\\f\\ L ^ m) 
The function u m = v m — u Sm v solves the problem, and we have: 

||Um||ff m < \\Vm\\H m + H^lk™ < \\f\\tf{n m ) + ||^V|| Am 
< ||/||La(n m ) + IMU' m < C||/|| L 2 (nm) 

3 A two-sided Dirichlet- Neumann domain de- 
composition algorithm for the Laplace op- 
erator 

Proposition 6 Let f G L 2 (Q) and u G Hq(Q) be the unique solution of 
Au = f. For m — 1,2 let f m = f\n m and g m G ifo(f2 m ) ^ e unique 
solution of Ag m = f m . Let i] m = g m \ T 

X = u\ r (5i + S 2 )X = -(Sir)! + S 2 rj 2 ) 
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proof: the direct implication is stating continuity of u and its normal deriva- 
tives through T. The converse implication states that taking u\Q m = 
solves the global problem. 

We use the same notation as in the previous proposition to state: 

Theorem 4 LetO < 9 < 1 with \\Si 1 S 2 + S 2 1 Si\\ C (a) < ^ ■ Any sequence 
(X n ) n C A which fulfills 

A n+1 = ((1 - 9)Id - 6 -(S^S 2 + S^Si))\ n - °-{S^ + S^mS^ + S 2V2 ) 

do converge in A (with geometric rate 1 — y at least) and its limit is u\r- 

proof: theorem [3] states selfadjointness of Si S 2 + S 2 Si in A and theorem 
Estates coerciveness of Si 1 S 2 + S 2 Si so: 

H(i _ )id - ^(Si's, + s 2 1 s 1 )\\ ci A ) = 

|(((1 _ e)Id - l(Si l s 2 + S 2 l S 1 ))\, A) A | 
a II a IIa 
(((1 - 9)Id - l(Si l S 2 + S 2 l Si))\, A)a V ( ^n^^, 36 ^ i 

T <i-0--(i+co<i- y <i 

4 Tools for a Domain Decomposition algo- 
rithm for the Helmholtz equation 

In the sequel we shall assume connnectedness of the open sets Q m , m = 1,2. 
4.1 On the j operator 

We define the linear bounded operator j from A to A' as the composition of 
the bounded linear injections 

a c h!(t) c L 2 (r) c fH(r) c a' 

It will play a key role for Helmholtz equations. Its properties are summarized 
by: 
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Proposition 7 . 

1- j is a compact and one-to-one operator 

2- For any A G A, j(X) = j{X) 

3- Form = 1,2 

VA G A, a G A, (j(A), (t)aa' = (<r, S m 1 j(X))A m = 0, X) Lm 
VAgA, (j(A), A) a ,a' = ||A||| 2(r) 

4- f = 3 

5- For m = 1, 2 t/ie operator S~ x j is selfadjoint in A m 

6- For m = 1, 2 t/ie operator jS~} is selfadjoint in A' m 

i 

proof: Item 1 comes from the Rellich compactness of the injection from H§ (T) 
to L 2 (r). 

Item 3 comes from: 

(a,S m 'j(X)) Am = J nm Vu^u s -j(X) = I r cr^^ = f r a\ = (*,\) LH n m) 
Item 4 follows item 3 because (j(A), ct)aa' = f r aX = (j(cr), A)aa' 
Item 5 follows items 3 and 4 because 

(v,S m 1 ](\)) Am = 0'(A),cr) AA / = (j(a),X) AAI = {S^j{a),X) Am 
Item 6 follows item 5 because 

(j^V^)a^ = (^J^V.C^a™ = {S^-^S^-jS^-v) Am = {n,jS^-v) K 

4.2 The spectrum of the local Helmholtz problems 

This paragraph is devoted to the study of the operators 

m = 1, 2 (S m + i 7 j) : A — ► A' 

related to the Laplace equation, and to the like (S^ +iyj) operators related 
to the Helmholtz equation. Let 7 denote a real number. 

Proposition 8 . 

1- VA G A ((S m + ijj)X, A)aa' = Pllim + ^7ll^lU 2 (r) 
z 7 j) /ias a bonded inverse. 

3-VAgA + iyj)X = (S m - z 7 j)A and T£X = T m 7 (A) 
proof: 

1- is straightforward applying proposition [7] and notation [2] 
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2- item 1 shows that Ker(S m + ijj) = {0} and Im(S m + ijj) is closed in 
A'. It remains to show that Im(S m + iyj) is everywhere dense in A'. By 
propositions [3] and [7] 

VA G A, ((S m + i 7 j)A, i])aa' = =^ 

VA, ((S m + i7j)?7, A)aa' = ==>> (S™ + i~/j)r) = ==>> rj = 

because (i> m + ryj) is one to one. 

3- is straightforward. 

Proposition 9 For m = 1, 2 and any z/ G A' there exists a unique u G if m 

9n r 
Am = m i2 m and — h ijjp u = v 

In fact u = u^ Sm+ilj " > ~ lp and 

\\u\\ Hm < C\\v\\ A > 

proof: Uniqueness is straightforward by the Green formula. For existence we 
apply the previous proposition to get A G A such that (S m + ijj)X = v, and 
check that u = u x solves the problem. 

The following remark will be crucial to prove convergence of domain decom- 
position algorithms for the Helmholtz equation: 

Remark 5 with the notation of the preceeding proposition, if v G L 2 (T) then 

- — G L 2 (T) and ||- — \\ L 2( r) < C||^|| L 2 (r) 
proof: -^- = v- i-fjp r m u G L 2 (T) + A C L 2 (T) and 

du . I. r i 

||^— \\l^(t) < IMk 2 (r) + |7ll|p m w|| L 2 (r) < \W\\l^[t} + bilMk„ 

< IMU 2 (r) +C|7|||^I|a' < C|MU 2 (r) 
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Proposition 10 For m = 1,2 and any f G L 2 (Q m ) there exists a unique 
u G H m such that 

du 

Au = f in VL m and h iljp V m u = 

on m 

and we have the estimate 

\\u\\ Hm < C||/|| L 2 (nm) 

proof: let v G .f/g(fi m ) solve Av = f. Let w = u — v. The function w 
fulfills Aw = and #^ + ^ ?>!> = Remark 1 shows that ^ G A'. 
The result follows from the previous proposition and the estimate follows 
estimates in remark H] and proposition [9j 



Remark 6 Using the notations of proposition [77J we have, as in remark 0, 
-Jr*- G L 2 (T) and the estimate: 

II ~— IU 2 (r) < C||/|| L2(Qm) 

(J I i m 

We can now proceed to compute the eigenfrequencies of the local Helmholtz 
problems involved in the Domain Decomposition algorithm. For that sake, 
we will use the following 

Notation 3 For m = 1,2 we denote by 

Di-.L 2 (n m )-^L 2 (n m ) 

/— ►« 

where u is given by proposition [77J 
This map has the following properties: 

Proposition 11 . 

1- is a compact operator in L 2 (Q m ) 

2- the adjoint map of for the L 2 (Q m ) scalar product is .D~ 7 

3- we have D^g = D^g and D^D^J = D^Dlf 

4- ImDl C H m 

5- o^Dlf G L 2 (T) 
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proof: 

1- translates Rellich compactness of the imbedding of H m in L 2 (Q m ) 

2- If u = D^f and v = D^g then Green formula gives 

/ ugdx— [ vfdx = [ J^- — udo— [ ^ U vda = —27 [ uvda+i^ [ uvda = 
Jn m Jn m Jr on m Jr on m Jr Jr 

3- is straightforward 

4- follows proposition [10] 

5- follows remark [6] 

We collect the spectral properties of in the following 

Proposition 12 We denote by a the spectrum of an operator and by a p the 
set of its eigenvalues. We have for m = 1,2: 

1- a(Dl) = {0} U a p (Dl) 

2- fie a(Dl) ^ /I G cr(-D m 7 ) 

3- For any f G L 2 (Vt m ), if u = D^f then 

(D^f, /)i2(n ro ) = - / I Vu\ 2 dx + 17 / \p v m u\ 2 da 

4- there exists a constant c such that If ' /i G a(D^ 7 ) ! // ^ 0, then 

Re /i < 0,7/m/i G R + , \Im ji\ < c\j\ \Re fx\ 

5- If 7 ^ i/ien 

) n R = {0} 

proof: 

1- follows compactness of D^ m asserted in proposition [TTJ 

2- is obvious by taking the complex conjugate of the eigenfunction associated 
with \i. 

3- 

( D mf\ /)i 2 (fim) = / ufdx = / uAudx = 



— f \S7u\ 2 dx+ [ p T m u do = — \ \S7u\ 2 dx + 27 / \p T m u\ 2 da 

Jn m Jr on m i r »/n ro 

4- if \i is an eigenvalue of with associated eigenfunction /, and u = D^f, 
then 

u / |f| 2 <ia; = — / \u\ 2 dx + w / Ip^upda 



18 



and the result follows with c the constant of continuity of the trace operator 
p r m from H m to L 2 (T). 

5- If \i 7^ is a real eigenvalue of with associated eigenfunction /, then the 
previous formula shows that p v m u = and, because u £ Ira(D^), = 

—ijp^u = 0. On the other hand, the equality D^f = pf translates to 
Am = -u. Because the Laplace operator is hyperbolic in the direction n m , 
and both data on T are zero, this implies u = on a neighbourhood of T. 
Solutions of elliptic equations being analytic, and fl m being connected, this 
implies u = on Q m . Then / = 0, which contradicts the assumption on / 
as an eigenfunction. 

Theorem 5 Let k 2 £ R, k 2 ^ 0. Let 7 £ i?, 7 ^ 0. For m = 1, 2: 

1- /or any / in L 2 (Q m ) and v £ A' , there exits a unique u £ H m such that: 

. , 2 „ du . . r 
Au+k u = f i1Jp m u = v 

OTlfYi 

2- we have the estimate 

\\u\\ Hm < C{\\f\\LHa m ) + \W\\a>) 

3- if v £ L 2 {T) then £ L 2 {T) and we have the estimate: 

du 

II \\v{v) < c{\\f\\ L 2 {Um) + IWWlhd) 

proof: 

1- Let A = (S m — i r yj)~ 1 i i and let v — u — u x . Then: 

a 7 2 r , du . . r 

Am + k u = j and — ffjp m u = v -<=>- 

on m 

dv 

Av = f — k 2 u x — k 2 v and — i'lip^Lv = -<=>- 

dn m 

v = DUf ~ k 2 u x - k 2 v) ^ Dlv + k~ 2 v = k~ 2 Dl(f - k 2 u x ) 

The previous proposition [12] shows that — k~ 2 £ a{D'J n ) so this problem is 
well-posed 

2- we have the estimate: 

\\v\\ LH n) < C(k- 2 \\Dlf\\ L2{n) + HD^II^n)) < C(Ar 2 ||/ || i2(n) + ||M A || L2(n) ) 
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so 

IMIl» ( o) < c'(k- 2 \\f\\ LHn) + \\u x \\ LHn) ) < C'(k- 2 \\f\\ LHn) + \\u x \\ H J 
< C'(k- 2 \\f\\ LHn) + ||A|| A ) < C'(k- 2 \\f\\ L 2 m + \\u\\ A .) 

We write 

a i- ,2 ® u • ■ r 

Au = f - k u - — - %!3p m u = v 

on m \ T 

and the H m estimate follows proposition [9] and proposition [KB 
3- If v G L 2 (d) we write again 

Au = f - k u - — - %!3p m u = v 
un m | r 

and the estimate follows remark [5] and remark [6j 



4.3 Despres operators and the energy fluxes 

We now define the building blocks of the intertwinning operator on the ficti- 
tious boundary T: the Despres operators. 

Definition 1 : Let k G R, 7 7^ 0, 7 G -R. For any z> G A' ; Zet u G -£f m fre 
i/ie unique solution, given by Theorem of the following equation on Q m , 
m = 1,2: 

2 ^ r 
Am + k u = 0, i^jp u = v on F 

dn m 

Let P^ be the linear bounded operator in A' defined by: 

67 5m , • ■ r r 

P rn V = a + ? 7JP m M O"- r 

OTl m 

Remark 7 Boundedness of P^ follows Remark^ and Proposition^ 
Proposition 13 : We obviously have P^ l P~' y = P'^P^ = Id^ 
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Notation 4 Let A 1 denote the linear bounded operator in A' x A' given by: 

Remark 8 : The inverse of A 1 in A' x A' is the bounded linear operator 
given by: 

(*>- - ( J r , "f ) 

In order to use conservation of energy, and to gain compactness, we use 
Theorem O to introduce: 

Notation 5 For 7 7^ and k G R, let the bounded operator in L 2 (T) denoted 
by P^ be the restriction of P^ to L 2 (T). Let the bounded operator in L 2 (T) x 
L 2 (r) denoted by A 1 be the restriction of A 1 to L 2 (T) x L 2 (T). 

Conservation of energy fluxes through V reads: 

Proposition 14 Let 7 7^ and m = 1, 2. 

(i) P^ is an isometry in L 2 (T): 

Vi/ e L 2 (T), |M| L 2 (r) = \\Ply\\ L 2 {T) 

(ii) A 1 is an isometry in L 2 (T) x L 2 (T): 

V(i/,77) G L 2 (T) x L 2 (r), \\^,r])\\ L 2 {rxL2(r)) = \\A^{u,r))\\ L z {T)xLHT) 
proof: For v G L 2 (T) let u G H m solve by Theorem the following equation: 

Au + k 2 u = in Cl m , — ijjp^Lu = v T 



Multiplying by u the equation fulfilled by u and integrating on Q m gives 

/ \S7u\ 2 dx — k 2 I \u\ 2 dx = f ^ U uda 
Taking the imaginary part gives 

f du _ 
Im / — — uda = 

Jr dn m 
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The result follows integration of the following identity on T: 

< du 2 du . . , 2 _ du 



An important consequence of this property will be crucial in the next section: 



Corollary 6 : For m = 1, 2 and 7 / 0, 

(%) P^ is a normal operator in L 2 (T) 

(ii) A 1 is a normal operator in L 2 (T) x L 2 (T) 



4.4 Spectral properties of the Despres operators 

In the preceeding section we proved that the Despres operator P£, ( m = 1,2 
and 7 7^ 0) is a bijective isometry in L 2 (T), and consequently a normal 
operator in L 2 (T). It follows that its spectrum is a subset of the unit circle 
in the complex plane. We now investigate this spectrum more accurately. 

Definition 2 Let 7 / and m = 1,2. Let be the operator in L 2 (Y) 
given by: 

WGL 2 (r), Clv = 3 p T m u 
where u G A is the solution given by Theorem [3] of the equation: 

(A + k 2 )u = in f2 m , iip v m u = v on V 

on m 

Compactness of the injection j from A to L 2 (7) gives: 
Proposition 15 : 

(1) For m = 1, 2 and 7 ^ ; P£ = I + 2i-fC^ 

(ii) is a normal and compact operator in L 2 ^) 

Notation 6 Let 

(i) denote the spectrum of P^ 

(ii) denote the sequence of eigenvalues of the Laplace operator on Q m 
with Dirichlet boundary condition on dQ m , i.e. 

—k 2 ^ S^ ir if and only if the following problem is well posed: 

(A + k 2 )u = f eL 2 (n m ), ueHl(Sl m ) 
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We have: 



Proposition 16 : Let 7 ^ 0, and m = 1, 2. 

(i) 1 belongs to 

^ 1 is an eigenvalue of P^ if and only if —k 2 £ S^ ir 

(i) Because is a compact normal operator (Proposition [13]) in L 2 (fl m ), its 
spectrum is a sequence of eigenvalues and its limit zero. This implies that 1 
is in the closure of which is closed. 

(ii) If —k 2 £ S^ ir , then there exist an eigenfunction ip^ of the Laplace oper- 
ator such that 

A<y9fc + k 2 Lfk = and Pk\r = ^ 

Let Uk = l^p- We have ^ ^ or else solving an elliptic equation, 
condition p^ = and connectedness of Q m would imply = on Q m , 
which contradicts the fact that pu is an eigenfunction. Obviously Uj. is an 
eigenvector of P^ for the eigenvalue 1. 

Conversely, if 1 is an eigenvalue of P^, with eigenvector v 7^ 0, then there 
exist u in H m such that 

Au + k 2 u = and v = — — — i^jp v m u = — — + i^jp v m u 

on m |p on m | r 

This implies p v m u = 0, which added to the fact that u £ H m implies u £ 
Hl{VL m ). So m is an eigenfunction of the Laplace operator for the eigenvalue 
—k 2 , provided u is not identically zero. And this is ruled out because = 
v ± 

The following Theorem gives a complete spectral description of the Despres 
operators: 

Theorem 6 : Assume that —k 2 ^ For 7 7^ and m = 1,2: 

^ £^ = {1} U (e l(7 ™) n (zN where (cr^) ng 7v is a sequence of numbers with 
o"^ £ R, 7^ ; and — > when n — ► 00. 

(ii) For each n & N, e %<T ™ is an eigenvalue of P^, with finite multiplicity. 
(Hi) L 2 (T) is the Hilbert direct sum of the eigenspaces associated with the 
eigenvalues e* ™ 
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Proof: is a normal and compact operator in L 2 (T) (proposition [To"]) . 
Its kernel is trivial (proposition [TBT) . The diagonalization theorem gives a 
sequence (A^) ne jv, (A^ 7^ 0), for its eigenvalues (they have finite multiplic- 
ity), with limit zero, and L 2 (T) is the Hilbert direct sum of the associated 
eigenspaces. The diagonalization of P^ = I + li^C^ follows, with eigen- 
values (1 + 2ryA^) ng jv, (1 + 2ryA^ 7^ 1). Proposition [T41 implies that these 
eigenvalues have modulus one: we set 1 + 2ryA^ = e*°™. 

4.5 Spectral Properties of PlP^ and P^P\ 

Properties of the intertwinning operator A 1 rely heavily on the spectral prop- 
erties of P\P^ and P^P\ that we investigate now. 

We first list obvious properties which follow from the previous section: 
Proposition 17 : 

(i) P\P~2 and P^P\ are isometric bijections in L 2 {^). 

(ii) P1P2 and P2P1 are normal operators in £ 2 ( 7 ). 

(Hi) P1P2 — / and P2P1 — / are compact operators in £ 2 ( 7 ). 

Proof: 

(i) follows the fact that P± and are isometric bijections in L 2 (T). 

(ii) follows (i) 

(iii) follows proposition [T5l through: 

p?pj = (/ + 2i 1 Cl){I + 2i 7 C 2 7 ) = / + 2vyCl + 2* 7 C 2 7 - 4 7 2 C 1 7 C 2 7 



An important property that we sail need is the spectral status of 1: 

Proposition 18 : 1 is not an eigenvalue of P1P2 or P^P\ in L 2 (T). 

Proof: By symmetry, it is enough to prove it for P^P^- Let v G L 2 (T) be such 
that P1P2V = v. This translates to the existence of U\ G H\ and U2 G H2 
satisfying: 

Au 2 + k 2 u 2 = in Q 2 ; 7— - — iku 2 = v, — - + iku 2 — Plv on T 

on? on? 
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v , . . r , du\ ., „, dui ., 
Aui + k u\ — m — %ku\ = P 2 v, — V iku\ — v on 

This implies that on V these functions fulfill: 



dui du 2 .. _ dux .. du 2 .. 

ikui - iku 2 = 0; h ikui - h iku 2 

oni on 2 on\ on 2 

Adding and substracting gives: 

du\ du 2 



u\ = —u 2 and 



drii dn 2 



We define u on Q as u\q ± = u\ and u\q 2 = —u 2 . It solves the Helmholtz 
equation on Q± and Q 2 , its has no jump accross T, neither has its normal 
derivative. So it solves Helmholtz equation on Q. Moreover u\qq = 0. As- 
sumption (A) gives u = 0, so U\ = and u 2 = 0; and v — follows. 



Proposition 19 : 

(i) The spectrum of PlP 2 in L 2 (T) is {1} U (e"~" 2 ) ng Ar ; where (r" 2 ) n6 Ar is an 
infinite sequence of real numbers, t™ 2 ^ 0, and r™ 2 — > when n — ► oo. 

(ii) (e lT ^) n( zN is the set of eigenvalues of P^P 2 - They have finite multiplicity. 
If we denote by E™ 2 the eigenspace associated with e""" 2 , then L 2 (T) is the 
Hilbert direct sum of the subspaces {E™ 2 )n&N 

(Hi) P 2 P\ has the same properties, and we set the obvious notations: (e lT ^) neN 
for eigenvalues and (E 2 i) n eN for eigenspaces. 

Proof: The operator P^P 2 — / is a normal compact operator (proposition 
ITTj) . So by the diagonalization theorem its spectrum is the union of {0} and 
an infinite sequence of eigenvalues with finite multiplicity (t n ) neN , {t n ^ 0). 
Zero is not an eigenvalue of P\P 2 — I (proposition [T8|) . So the whole set of 
eigenvalues is {t n ) n& N- If E\ 2 denotes the eigenspace associated with t n , then 
L 2 (T) is the Hilbert direct sum of (E™ 2 ) neN . By proposition [T71 we know that 
P\P 2 is an isometry in L 2 (T), so |1 + t n \ = 1, and we write it: t n = e lT ^. 
The theorem translates proven properties of t n into properties of r™ 2 . 

In order to study the relationship between (rf 2 ) ng Ar and (r 2l ) n€N , and be- 
tween (E™ 2 ) neN and (E 21 ) neN , we prove the following lemmi: 
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Lemma 1 Let 7 7^ 0. For m = 1, 2 and any f G L 2 (T): PmV = P^v 
Proof: Let m = 1 or 2. By definition of P^ there exists n m G H m with: 

Au m + *V = 0; — - % 13 p v m u m = u; j— + i 1 3 p v m u m = P^v 
on m | r an m | r 

Taking the complex conjugate of these equalities gives n™ G -ff m such that 

A?r + /cV^ = 0; - — + ijjpLvT 1 = V; - — - i^jpLW 1 = P&u 
on m \r on m \T 

which by definition of P^ writes 



P-^V = Plv 

Lemma 2 Let 7 7^ 

(i) If X is an eigenvalue of P^P 2 (resp P 2 P±) for the eigenvector v then it 
is an eigenvalue ofP 2 Pi (resp PlP 2 ) with associated eigenvector V 

(ii) For all n G iV , r™ 2 = t 21 {mod 2tt); we denote it by r n 

(Hi) If we denote by C the set of complex conjugates of distributions in a set 
C , then, for any n G N, 

Wi = E? 2 and M = E& 

Proof: 

(i) If P^P 2 v = \v then PlP 2 v = Xu which by lemma Q] writes P 1 ~ 7 P 2 ~ 7 i 7 — 
Xn, which implies, by proposition [T3l P 2 Plv = =u = Xu because = = A by 
proposition [T9J 

(ii) follows (i) and a renumbering. 

(iii) follows (i) because it gives: 



#21 C E^ 2 and C 

but then 



E&CE&CE*! and E? 2 C E& C E? 2 

which gives 

Wi = K and Ej 2 = E n 2X 
Lemma 3 ; Let 7 7^ 0. For any n e N, P^E^ = E? 2 and P^E^ = 
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Proof: Let c^Oe then P?Pjv = e iTn u so P 2 7 P 1 7 P 2 7 z/ = e iT "P> which 
proves that P^v £ P 21 (-^^ 7^ because P2 7 is bijective (proposition [TBI)). 
This writes P^E^ C P 21 . Proposition [T9l lemma [2] and invertibility of P^ 
(proposition [13]) give dimE^ = dimE™ 2 = dimP^E^ so P 2 7 P™ 2 = P 21 . 

The following algebraic property and its consequences on the eigenprojectors 
(next theorem) are a key for understanding the geometric properties of the 
intertwinning operator: 

Lemma 4 : Let 7 7^ 0. For any /i ^ {1} U (e lTn ) ne N-' 

(pm - /ipr 1 ^ 7 = nwn - M-0 -1 

(P 7 P 7 - /x/^P 7 = P 7 (P 7 P 7 - /x/)- 1 
Proof: we prove the second assertion, using resolvant identity: 

(P 7 P 7 - /i/)-^ 7 = (P 7 P 7 - rt^PSWPS - /i/)(P x 7 P 2 7 - ^PT 1 = 

[[/+Ai(W-A0"W-A*(^-A^- 1 ^](W-M0- 1 = P 2 7 (P7P 2 7 -///)- 

Lemma 5 Le£ 7 ^ 0. Por any n G iV, i/n" 2 (resp. II 21y ) denotes the spectral 
projector of the operator P^P^ (resp. P%P\) on £/ie eigenspace P™ 2 (resp. 
E21) then we have 

p x 7 n^ = n? 2 p x 7 and p 2 X 2 = n^p 2 7 

Proof: By symmetry it is enough to prove the first formula. Let C n denote a 
positively oriented curve in the complex plane, which winds one time around 
the eigenvalue e %Tn , and none around any other eigenvalue, then the Dunford 
integral representation formula gives: 

n? 2 = ^ j Cn ( p i p 2 - mO" 1 ^ ^d n£ = ^ / Cn (W - »iY l dii 

The previous lemma gives the following: 

PIK2 = ^f c PKPlPl-vir'dv = ^f c (P 2 7 Pi 7 -///)- 1 P 2 7 rf/i = n^p 2 7 
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4.6 Spectral Properties of A 1 

We recall that 

A 



_ -P 7 



and that this operator is a bijective isometry in L 2 (Y) x L 2 (Y) 
Theorem 7 : Let 7 7^ 0. 

(%) // A ^ {±1} U (±e*^ L ) ne 7v ^en A belongs to the resolvant set of A 7 and 

M7 _ in -i_f A(P7P 2 7 -A 2 /)- 1 -(P1P2 - \ 2 I)~ l Pl \ 
{A M} ~ y -(P2P1 - A 2 f) _1 P 2 7 A(P 2 7 P 1 7 -A 2 /)- 1 J 

(m,) For any n £ N, ±e l ~%~ is an eigenvalue of A 1 with associated eigenspace: 

F± = {(jjl , Te"^P 2 7 /i); /i G E&] (2) 
and associated eigenprojector : 

fiiij {±1} belong to the spectrum of A 1 and are not eigenvalues of A 1 

(iv) (F^) neN .± is an orthogonal family of sub spaces and we have the Hilbert 
decomposition 

L 2 (r) x L 2 (r) = (© °°f+) © (®-p-) 

(v) The following series are strongly convergent in C(L 2 (T) x L 2 (T)): 





00 00 

^ = £e^ + -£e^P n - 



(i) let A^{±l}U(±e^) neJV 



Proof: 
(i) let 

(A 1 - XI) is injective: let (ip,ip) e L 2 (T) x L 2 (T) be such that 

' ip \ _ ( \ 



(A 1 - XI) 
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this writes 

PH) + Ay? = and P^tp + = 

which implies 

P 2 7 P 1 7 ^-AV = and P%P%ip - X 2 (p = 

This implies y? = ip = by proposition [19] and lemma [2j 
(A 7 - XI) is surjective: For any (£,77) e L 2 (T) x L 2 (T) let: 

V = (PM - A 2 /)- 1 (AC - P7??) and ^ = (A 7 A 7 - A 2 /)" 1 (Ar / - P 7 £) 
We have, by lemma H] 

P 2 > + A^ = P 2 7 (P 1 7 P 2 7 - X 2 I)-\H - p iV) + A(P 2 7 Pi 7 - A 2 J)- 1 (At; - P 7 £) 

= (P 7 P 7 - A 2 I)- 1 (AP 2 7 £ - P 7 P 7 r/) + A(P 7 P 7 - A 2 !)" 1 ^ - P 7 
= (P 7 P 7 - A 2 /)- 1 (-P 2 7 P 1 7 7 7 +A 2 t 7 ) = -77 

Similarly 

P7^ + A^ = -e 

These two equalities write 

So surjectivity is proven. These expressions for ((p, ijj) give the formula for 
the resolvent of A. 

(ii) By definition of P™ 2 we have for any \i G P™ 2 : 

^(/i , T e-^P 7 /i) = (±e-^P 7 P 7 /i , -P 7 /i) = 

(±e^/i , -P 2 » = ±e i2 *(/x , T e-^P 2 7 /i) 

Because P™ 2 7^ {0}, this proves that ±e*~ is an eigenvalue of A 1 . This 
proves moreover that is a subset of the eigenspace of A 1 associated with 
the eigenvalue ±e i_ 2\ 

On the other hand, if (£, 77) is an eigenvector of A 1 for the eigenvalue ±e l2 ? 
then 

-P 7 77 = ±e^£ and - P 2 7 £ = ie^ry 
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This implies 

P?P?t = T* i! *P?T, = e* r »t so £ e 

and 

77 = Te^PK 
This completes the caracterisation of the eigenspace. 

We compute now the eigenprojector: for this sake, we make a choice of a 
branch for yfz. We take a positively oriented curve in the complex plane 
which winds one time around ±e*^ and not around ~^e % ~^ nor does it wind 
around any ±e l_ 2~ for n' ^ n. Let D n be the image of C„ by the function 
z — > z 2 . L> n winds one time around e* rn and does not wind around e lTn ' for 
n' 7^ n. Let Z)^ wind one time around e lTn , lying in the interior set delimited 
by D n . 

The eigenprojector is given by the Dunford formula: 



p± = — I (a 1 - xn^dx 

2in Jcr ' 



Using the representation formula given by (i) for (A 1 — A/) -1 leads to compute 

integrals of two different types: 

For the first type it is straightforward and gives: 



For the second type, we first use the resolvant identity to have: 

^tL^- XI ^x ^L}^-^= 

I I ( p i p 2 ~ A/)~ 1 (P 1 7 P 2 7 - /i/)" 1 ^ = 
2m2mJ Dn JDC 12 ; v 1 2 P ; 2^X 

2m2mJ Dn y 1 2 ; \J D> n X - pi J 2y/X 

- — — f ( W - /i/)^ f / ) dp = 

2m2mJDC 12 ' \J D n 2VX(X - pi) J * 
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2m V 12 ^ > 2JJi 

We compute now the second type of integral using this equality, properties 
of Il™ 2 , and lemma [5] to have: 

- xmdX = - L {p ^ - A/r1 ^) p " = 

-(—[ (p?pZ - xir 1 ^^) iq = - [ , . 1 n? 2 -^=W 

\2mJ Dr } 1 2 ; 12 2v / Ay' \2m J Dn {e^-XI) 12 2V\) 1 

(iii) ±1 are limits of the sequence of eigenvalues (±e^ IL ) ne Ar so they belong 
to the spectrum of A. These values are not eigenvalues, or else 1 is an 
eigenvalue of P^Px and PiP 2 , which is ruled out by proposition [T8l 

(iv) and (v) Assertions (i), (ii) and (iii) prove that the spectrum of A 1 is 
{±1} U (±e*^") ng 7v- Normality of A 1 (corollary [6]) implies orthogonality of 
the family (F^) n ^N, and gives the decomposition of / and A 1 as series of the 
eigenprojectors (P„ ) 

Remark 9 : Notice that the expression of in (ii) of the previous propo- 
sition is symmetric: in fact we have 

{(ji , T e-^P^); fi G E? 2 } = {(zpr^Pfy , //); // G E%} 

this is because PxE 21 = E\\ following lemma [3], so if we set ji = =Fe~ l ^ L P 1 y / i ' 
it ensures fi G E\\ if A*' G E 21 . Moreover, by definition of \J we have: 

Te^P^fi = e-^PjP^' = fi' 
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5 Domain Decomposition algorithm for the 
Helmholtz equation 



5.1 The Domain Decomposition framework for the Helm- 
holtz equation 

Proposition 20 : Let k £ R, k ^ 0. For any f £ L 2 (Vt), let u £ H^(Q) be 
the unique solution of the Helmholtz equation 

Au + k 2 u = f 

Let 7 £ R, 7 ^ 0. For m = 1,2, let v m £ H m solve the equation 

dv 

Av m + k 2 v m = /|n m and — — - i~fjp r v m = on T 

on m 

Let 

v m = 2i-fjp T v m and r] = (Pih> 2 , P^i) 
Then the equation in L 2 (T) x L 2 (T) 

(A 1 - Id)n = r] 

has a unique solution: 

it = (tti, tt 2 ) with vr m = h iljp T u - 2i>yjp r v m 

on m 

proof: First notice that the assumption / £ L 2 (Q) and assumption (A) on k 
imply u £ L 2 (Q), so Am £ L 2 (Q). Regularity of dQ enables the use of classical 
regularity results ([Ag]) for solutions of elliptic boundary problems to have 
u £ H 2 (VL), hence ^ £ H*(T) C L 2 (T). This proves that n m £ L 2 (r). 
Let w m = u\ nm - v m . Then 

du . . r r.- • r dw m . . r 

TTm = ~ 1" 11JP U - 2%'jjp V m = h VJJ P Wm 

on m dn m 

Because the function w m fulfills 

2 dw m r du r 

Aw m + fc w m = 0, w m £ if m , ryjp w m = i-fjp u on r 

on m dn m 
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one has 

or/ du r \ 

OTl m 

so, if m' = 2, 1 for m = 1, 2 

= -^m(^ + «7Jp r M) = -Pl^m' + I'm') 

This writes 

71 = A y 7l — 1] 

Uniqueness follows from theorem [7] (iii). 

Proposition 21 : Let k G R, k ^ 0, 7 G R, 7 7^ 0. For / G L 2 (fi) and /or 
m = 1, 2, let v m & H m solve the equation 

du 

Av m + k 2 v m = f\n m and — — - ijjp r v m = on T 

Let 

v m = 2i'jjp r v m and 77 = (Ffi^, P^i) 
Let 7r = (71*1,, 7T2) G -^ 2 (r) x L 2 (T) solve the equation: 

(A 1 - Id)n = n 

and let u m G H m solve the equation 

du 

Au m + k 2 u m = f\n m and — ^ + i^jp v u m = Tr m + v m on T 

on m 

Then u given by u\a m = u m solve the Helmholtz equation 

Au + k 2 u = f ueH^(n) 

proof: By definition of u m and v m one has: 

A(u m f m ) -|- k {u m fnj) 0, u m v m G H m 

and 



+ hjP («m -V m )=TX m On T 



<9n TO 

33 



This implies through proposition [Tj3 

/ m \-i -7 d(u m — v m ) . . r . <9w m . . r 

(P^J 7T m = P m 7 7T m = — r/jp (tt m - U m ) = — Ijjp U r 

Because 7r = A 1 it — r] this implies (with m' = 2, 1 for m = 1, 2) 

du>ffi p ^^tyi 1 p 
«7jp U m = -7T m / - V m i = -- i^jp u m > 

on m dn m > 
Adding and substracting these equalities gives: 
du m du m > 



and p r u m = p T u m i 



These jump conditions through V imply that Am + k 2 u = f on Q, and u 
fulfills the Dirichlet boundary condition on dfl because u m G H m . 

Remark 10 Theorem^ shows that the problem 

(A - Id)n = rj 

is ill-posed for 7] £ L 2 {T) x L 2 {T). Proposition [20] shows that if t] has the spe- 
cific form given through the domain decomposition setting for the Helmholtz 
equation, the equation (A — Id)n = r\ do have a solution, (and this solution 
is unique by Theorem [7|). Proposition \2J] shows that this solution provides 
the solution of the Helmholtz equation. 



5.2 The domain decomposition ^-algorithm for Helmholtz 
equation 

Let / G L 2 (Q). Let u £ Hq fulfill the non-dissipating Helmholtz equation 
Au+k 2 u = f in Q. The classical algorithm used, (for dissipating cavities with 
Sommerf eld-like boundary condition), to solve by a domain decomposition 
technique the Helmholtz equation ([B],[BD],[D1],[D2],[CGJ]) writes, in 
the non-dissipating case that discussed here, as follows: for any tt° = (7r°, n®) 
given in L 2 {Y) x L 2 (T) let 

n P+i = 07T p + (x _ e)A^ix p - (1 - 9)r] 
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where r\ = {P^Vi, P^i) with v m = 2i r yjp T v m for v m G H m solving 

Av m + k 2 v m = f\n m and — — - i-fjp r v m = on T 

on m 

It is straightforward to translate the ^-algorithm in a PDE setting: one use 
theorem [5] to get the (unique) function G H m such that 

Aw p m + k 2 w p m = and n p m = — ^ + i 1 jp T w p m 

on m 

In these w p = (w p , w p ) variables the ^-algorithm becomes: for m = 1 and 2 
(resp. m! = 2 and 1) 

A< +1 + fc 2 < +1 = 0, < +1 G# m 

^-^P r < +1 = fi[^-<7i/<]-(l-fi)[^+<7^< f +^] on T 
cn m an m op m > 



For practical use in computing codes, one writes this algorithm in the u n = 
(w™, u%) variables with u v m = w v m + v m and gets: 

AuC 1 + k'u^ 1 = f lnm , u p + 1 eH m 

^ - ijjfiC 1 = 8[P^ ~ iljp r <] ~ (1 - e)[%&- + H3P T u p m ,] on Y 
on m on m on m ' 

5.3 Convergence results for the ^-algorithm 

Notice that if the sequence (7r p ) pe 7v has a limit 7r°° in L 2 (T) x L 2 (T) then 
continuity of A 1 gives: 

(A 1 - Id)ir°° = rj 

and 7r°° provides the solution of the Helmholtz equation on Q as stated in 
proposition [2Tl 

Alternatively, a way to solve the Helmholtz equation on Q through solving 
Helmholtz equations on Q m , (m = 1,2), is to notice that convergence of 
(7r p )p(zN in L 2 (T) x L 2 (T) implies convergence in A' x A', and theorem [5] shows 
that the sequence (w p ) p ^n = ( {w p , w p ) ) p ^n has a limit in H\ x H 2 , which 
implies convergence of the sequence (u p ) p< z^ = ( (u p , u p ) ) pe jv in H\ x H 2 . 
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Proposition I2T1 shows that its limit u°° = (u^°,u'^ > ) provides the solution of 
the Helmholtz equation, through. 

{u™,u™) = (u|ni,W|n 9 ) 

Here are the convergence results for the ^-algorithm. We begin with a nega- 
tive result: 

Proposition 22 If 9 = then the sequence (7i p ) pe ^ has no limit in L 2 (T) x 
L 2 (r) unless if its initial value fulfills {A 1 — Id)ir = rj. Written for the 
(u p ) p£ N sequence, this translates to u\ = U\n t andu® = u\^ 2 where u G Hq(Q) 
solves Am + k 2 u = f in Q 

proof: for 9 = one has ir p — n p ~ 1 = y4 7 (7r p_1 — 7r p ~ 2 ) and proposition 
gives 

Vp ||7T P - 7r P_1 || L 2( r)xL 2( r) = ||7T P_1 - 7T P_2 || L 2( r)xL 2 (r) 

This prevents convergence unless if tt 1 = 7r°, i.e. ir° fulfills 
i.e. unless v,® = and u 2 = u\n 2 by proposition [2T1 

Remark 11 If 9 = the sequence (u p ) p£ n may have a limit in Hi x H 2 
even if u\ ^ u\q x or u 2 ^ tt|n 2 . This is because convergence of (u p ) peN in 

Hi x H2 implies convergence of + 27 j p l u\, + i / yjp r u 2 ) in A' x A', i.e. 

convergence of + i^jp v u P i ) ^| + i"ijp v u^) in A' x A' ; ie. convergence of 
(7r p )p g Ar m A' x A' ; which do not contradict divergence in L 2 (T) x L 2 (r). 

We now turn to the main result: 

Theorem 8 : For f G L 2 (Vt) and k G R let u G i/"o(0) so/fee Am + k 2 u = f. 
Let 7 7^ 0, 7 G i?. Then for any < 9 < 1 : 

(i) the sequence (it p ) p& n given by the 9 -algorithm converge in L 2 (T) x L 2 (T) 
to n u = «, vr£) with: 

9u , • ■ r 1 « Sit , . . r 

7r| = h ryjp tt and vr 2 = h t'JJp u 

oni on 2 

(ii) the sequence (u p ) p£ n given by the 9-algorithm converge in Hi x H 2 to 
(w[n is «|na) 

(Hi) There is no uniform geometric rate of convergence. 
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Proof: 

(i) let vr n = (71?, vr^) with: 

„ 5m . . r 1 •» ® u ■ ■ r 

7Ti = 7: h«7JP w and 7r 2 = — h«7jp w 

ani an 2 

Nullity of jumps of w and its normal derivatives through T gives 

vr" = Air u + r] 

This implies 



n p -n u = eh 11 - 1 - 7T U 1 + (1 - 9)A[ 7 i p - 1 - vr"! 



We use eigenprojectors of A 1 given by theorem [7] and denote by: 

Completeness of the set of orthogonal eigenprojectors (P*) n ,± proved in 
theorem [7] gives: 

7r"-7r w = + 

n n 

Decomposition of L 2 (T) x L 2 (T) by eigenspaces of A 7 writes for successive 
terms of the algorithm sequence as follows: 

This implies: 

IK±IU 2 (7)xL2( 7) = [1 - 20(1 - 0)(1 =F cos y)] p ||^ j± |U2 (T)xL 2 (7) 



and orthogonality of the eigenprojectors writes: 

-JJl ' II V±llL 2 ( 7 )xL 2 ( 7 ) 



Hl! 2 ( 7 )xL 2 ( 7 ) = Et 1 - 2^(1 - 9)(1 T cos ^)] 2p ||< ± ||i 2(7)xL2( 

n,± 



Assumption < 6> < 1 and theorem [7J (with proposition [18] asserting r n 7^ 
mod 27r) imply 

< 1-20(1 -0)(1=F cos—) <1 

2 

and Lebesgue convergence theorem gives 



L 2 (r)xL 2 (r) t 

7T ► 7T 
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This proves assertion (i). 

Assertion (ii) is straightforward: convergence of tt p to ix u in L 2 (T) x L 2 (T) 
implies its convergence in A' x A' which implies convergence of the related 
sequence w p in Hi x H 2 , and accordingly convergence of u p to (w^, «|n 2 ) 

Assertion (iii) is obvious by taking initial data for tt p in the n— th eigenspace 
of A 1 and notice that r n — > (proposition [T9l) 
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